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ABSTRACT: Large classes of new physics theories predict the existence of new scalar states, 
commonly dubbed sgluons, lying in the adjoint representation of the QCD gauge group. 
Since these new fields are expected to decay into colored Standard Model particles, and 
in particular into one or two top quarks, these theories predict a possible enhancement of 
the hadroproduction rate associated with multitop final states. We therefore investigate 
multitop events produced at the Large Hadron Collider, running at a center-of-mass energy 
of 8 TeV, and employ those events to probe the possible existence of color adjoint scalar 
particles. We first construct a simplified effective field theory motivated by i?-symmetric 
supersymmetric models where sgluon fields decay dominantly into top quarks. We then use 
this model to analyze the sensitivity of the Large Hadron Collider in both a multilepton plus 
jets and a single lepton plus jets channel. After having based our event selection strategy 
on the possible presence of two, three and four top quarks in the final state, we find that 
sgluon-induced new physics contributions to multitop cross sections as low as 10-100 fb 
can be excluded at the 95% confidence level, assuming an integrated luminosity of 20 fb -1 . 
Equivalently, sgluon masses of about 500-700 GeV can be reached for several classes of 
benchmark scenarios. 
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1 Introduction 

The Standard Model of particle physics has passed for many years all experimental tests. 
Only the mechanism of electroweak symmetry breaking, which is currently being addressed 
by the general-purpose experiments ATLAS and CMS of the Large Hadron Collider (LHC), 
at CERN, has remained for a long time an unsolved question. The recent observation of 
a neutral bosonic particle compatible with a Standard-Model- like Higgs boson [1, 2], to- 
gether with the imminent measurement of its properties, would then represent an impressive 
success of this theoretical framework. However, the mass of a fundamental scalar field is 
drastically affected by quantum corrections, which leads to the conceptual question of its 
stabilization with respect to the Planck scale lying orders of magnitude away from the weak 
scale. Over the last decades, large classes of alternative theories have been proposed in order 
to extend the Standard Model and cure this issue. Among these, weak-scale supersymme- 
try, and in particular its minimal version, the Minimal Supersymmetric Standard Model 
(MSSM) [3, 4], is one of the most theoretically and experimentally studied option. It is 
known to solve this so-called "hierarchy problem" by introducing partners to the Standard 
Model degrees of freedom with opposite statistics. In addition, several other conceptual 
problems of the Standard Model are addressed, such as the unification of the gauge cou- 
pling strengths at high energies or the question of a viable dark matter candidate. 

Experimental searches for the supersymmetric partners of the Standard Model particles 
are therefore among the main topics investigated at the LHC. Up to now, both the ATLAS 
and CMS collaborations have mainly focused on the strong production channels, since they 
yield larger cross sections. As a result, no traces of squarks and gluinos have been detected 
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so far and the limits on the masses of the first-generation and second-generation squarks as 
well as those on the mass of the gluino are consequently pushed to a higher and higher scale 
[5-7]. Therefore, the experimental attention starts to shift towards third-generation squarks 
and electroweak production channels. However, all the current results may not be valid for 
more general supersymmetric scenarios. They are indeed either derived in the framework of 
the constrained version of the MSSM, where the hundreds of free parameters of the general 
model are reduced to a set of four parameters and one sign, or in the context of simplified 
models inspired by the same constrained MSSM. In contrast, there are a vast variety of 
non- minimal supersymmetric models which are valuable to be investigated at the LHC. In 
particular, some of the final state signatures predicted by these non-minimal models require 
dedicated phenomenological studies in order to be ready for the interpretation of the data 
in these theories. This is one of the scopes of the work presented in this paper. 

Motivated by these arguments, we address the production and decay of color-octet 
massive scalar particles, also dubbed sgluons, that are predicted by several non-minimal 
supersymmetric models. The most considered examples are N = 1/N = 2 hybrid supersym- 
metric theories [8-16] and i?-symmetric supersymmetric theories [17-19]. In these models, 
the vector supermultiplet of the MSSM associated to the QCD gauge group is supple- 
mented by an additional chiral supermultiplet lying in the adjoint representation of SU(3). 
This new supermultiplet contains, on the one hand, a two-component fermionic component 
which mixes with the usual gluino to form a four-component Dirac fermion. On the other 
hand, it also includes a color-octet complex scalar particle, i.e., a sgluon field. Let us note 
that color octet scalar particles also appear in vector-like confining theories [20-25] or in 
extra-dimensional models [26]. 

As all of these models are similar in the sense of the LHC signatures, we adopt the 
approach introduced by the LHC New Physics Working Group [27] and employ a simplified 
model describing a scalar octet field and its interactions with the Standard Model sector 
[28]. This has the major advantages to leave open the possibility of reinterpreting the 
results in the context of any of the original models (or even in the framework of another 
theory including a sgluon field) and also to avoid handling a complete model which demands 
a careful design of a theoretically motivated but not experimentally excluded benchmark 
scenario. Equivalently, rather than fixing hundreds of free parameters related to one of the 
new physics theories above-mentioned, we focus on a specific sector of the model, relevant 
for our study, that is described by a small number of couplings and masses. 

We therefore start by constructing a simplified model describing the dynamics of a 
scalar field lying in the adjoint representation of the QCD gauge group. Motivated by com- 
plete models where loop-induced operators imply sgluon fields singly coupled to Standard 
Model quarks and gluons, we also include the corresponding interactions in our effective 
theory. Subsequently, once produced through standard strong interactions, a sgluon can 
then decay either to a quark pair or to a gluon pair, the latter being in general dominant 
for a sgluon field with a rather low mass [10, 12]. If kinematically allowed, the same inter- 
actions ensure a possible decay to a pair of top quarks or to an associated pair comprised 
of one single top quark and a light jet. 

In this scenario, sgluon pair-production and decay at the LHC could lead to signatures 
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possibly containing two, three or even four top quarks, the last two final states being 
expected to be largely suppressed in the context of the Standard Model. In this paper, 
we perform a phenomenological analysis of the sensitivity of the LHC collider, running at 
a center-of-mass energy y/s of 8 TeV, to the possible observation of sgluon fields through 
multitop signatures. This work extends our previous contribution to the 2011 Les Houches 
workshop [28] by including a more accurate description of the Standard Model backgrounds 
and a more efficient search strategy. It is also complementary to recent ATLAS and CMS 
investigations where sgluon-induced four-jet signatures have been considered [29-32]. In 
these analyses, limits on the sgluon mass up to 2 TeV have been extracted. In our theoretical 
setup, we evade these bounds by imposing the sgluon field to dominantly decay to at least 
one top quark, as motivated by realistic supersymmetric models including a color-adjoint 
scalar particle. Consequently, lower sgluon masses can be possibly expected. 

This paper is organized as follows. In Section 2, we describe our proposal for a simpli- 
fied theory modeling sgluon pair production and decay at the LHC and define benchmark 
scenarios for our phenomenological studies. For each benchmark point, we present the 
sgluon decay table and the leading-order and next-to-leading total production cross sec- 
tions. Section 3 is dedicated to our technical setup for the Monte Carlo simulation of both 
signal and background events at y/s = 8 TeV. A particular emphasis is put on the fast 
detector simulation package that we have used and in the way objects are reconstructed. 
In Section 4, we present the details of our phenomenological analyses and the associated 
results. Our conclusions are given in Section 5. 



2 A simplified model for sgluon production and decays at the LHC 



In order to investigate sgluon production and decay at the LHC, we construct a simplified 
model describing sgluon interactions with the Standard Model fields. To this aim, we follow 
the approach of Ref. [27] and extend the Standard Model in a minimal way by supplementing 
to its particle content one massive real scalar field a lying in the adjoint representation of 
the QCD gauge group. Its kinetic and mass terms are standard and can be expressed in 
terms of the QCD covariant derivative taken in the adjoint representation, 



£kin = -Da^D^da - \mla a a a with 



D„a a = d„a a + g s f bc a G* a 



(2.1) 



where adjoint gauge indices are explicitly indicated. In the expressions above, we have 
introduced the strong coupling constant g s , the antisymmetric structure constants of SU(3) 
fb c a , the gluon field and the sgluon mass m a . 

These interactions do not include single sgluon couplings to Standard Model fields. 
However, the presence of additional particles in general induces, in complete theories, loop 
diagrams leading to effective operators describing the coupling of a single sgluon field to 
up- type (u) and down- type (d) quark pairs, as well as to gluon pairs. In our theoretical 
framework, those interactions are described by the effective Lagrangian 



£ e ff = <? a dT a 



a L d P L + a*P R 



d + a a uT a 



atPL + a*P R u + agda^G^GV c + h.c. . (2.2) 



be a. 
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The matrices T a and the tensor d a are respectively the fundamental representation ma- 
trices and the symmetric structure constants of SU(3), the operators Pl and Pr are the 
left-handed and right-handed chirality projectors acting on spin space and G^ u a is the gluon 
field strength tensor. We have also introduced the parameters a q and a,?, with q = u,d, 
to model the strengths of the interactions among left-handed and right-handed quarks q 
and a single sgluon, respectively, as well as one single parameter a g for the modeling of the 
interactions among two gluons and one sgluon. It is important to note that the interactions 
included in the Lagrangian of Eq. (2.2) open all the possible sgluon decays to Standard 
Model colored particles. 

Inspecting the two Lagrangians of Eq. (2.1) and Eq. (2.2), one observes that our sim- 
plified theory is described by one mass parameter, the sgluon mass m a , and the effective 
couplings of sgluons to colored partons described by a set of four complex 3x3 matrices 
in flavor space a^, af, and ct„ , together with one complex number a g . Assuming 0(1) 
effective interactions, sgluon masses up to about 2 TeV are already excluded by dijet res- 
onance searches [29]. To evade this constraint, we choose scenarios where the sgluon field 
dominantly decays into final states containing at least one top quark or where its couplings 
to a pair of light quarks or to a pair of gluons are reduced. This choice is motivated by 
i?-symmetric supersymmetric theories and by N = 1/N = 2 hybrid supersymmetric models. 
In these models, interactions among a single sgluon and a pair of quarks are driven by loops 
of squarks and gluinos. The computation of those loops lead to effective couplings of the 
order of m q /M susy , M susy being the order of magnitude of the squark and gluino masses 
and m q the mass of the heaviest of the external quarks. These interactions are therefore 
non-suppressed only if at least one of the external quarks is a top quark. Furthermore, 
the coupling strength a g is related to a dimension-five operator and is thus expected to be 
suppressed too. One finds, from the explicit calculation of the relevant loop diagrams, that 
a g ~ 1/M susy [10]. 

We therefore consider two series of benchmark scenarios, refered to as scenarios of class 
I and class II in the following. For the first set of scenarios, sgluon particles are allowed to 
decay, in a universal way, to any associated pair of up-type quarks containing at least one 
top quark. We subsequently fix 

{a L u fi = (a*)\ = {a L u )\ = (a*)% = 0.3 , (2.3) 

for z = 1,2,3 and impose that all other interactions among quarks and sgluons are vanishing. 
We have in addition assumed that left-handed and right-handed quarks contribute in the 
same way to the interaction strengths, i.e., the left-handed and right-handed couplings are 
taken equal. Concerning the parameter a g related to the strength of the coupling among 
sgluons and gluons, we choose the value 

a g = 1.5 x 10~ 4 GeV" 1 . (2.4) 

Both Eq. (2.3) and Eq. (2.4) correspond to a scale M susy of 0(2 TeV), in agreement with 
the current experimental results related to the direct searches for squarks and gluinos at 
the LHC [6, 7]. 
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Parameters 


Scenarios of type I Scenarios of type II 


a g 
M 3 3 
(a u ) 3 i = (a«) 1 3 
{a u ) 3 2 = {a u ) 2 3 
m„ 
m t 


l.o x 1U Uev l.o x 11) (jev 
0.3 0.3 
0.3 
0.3 
[200-1000] GeV [400-1000] GeV 
172 GeV 172 GeV 



Table 1. Non-zero input parameters for benchmark scenarios of class I (second column) and II 
(third column). All the Standard Model parameters but the mass of the top quark (indicated in 
the Table) follow the conventions of Ref. [33] . 

In scenarios of class II, we focus exclusively on sgluon-induced LHC signatures with 
four top quarks and therefore forbid flavor violation in the sgluon decays. The only non- 
vanishing effective interactions are thus driven by the parameters 

(a£) 3 3 = (a,f ) 3 3 = 0.3 and a g = 1.5 x 10" 4 GeV" 1 . (2.5) 

The values of the non-zero parameters of the Lagrangians included in Eq. (2.1) and 
Eq. (2.2) are summarized, for the two series of benchmark points, in Table 1. For all the 
Standard Model parameters, but the top quark mass mt which is chosen equal to 172 GeV, 
we follow the conventions of Ref. [33]. Moreover, the sgluon mass is kept free and allowed 
to vary in the range m a G [200 - 1000] GeV and m a G [400 - 1000] GeV for class I and 
class II scenarios, respectively. 

A key element in the multitop analysis of sgluon production and decay at the LHC lies 
in the sgluon branching fraction to final states containing one or two top quarks. To study 
the evolution of these branching ratios with the sgluon mass, the Lagrangians of Eq. (2.1) 
and Eq. (2.2) have been implemented in the FEYNRULES package [33-37] and the model 
has been exported to the UFO format [38]. We have subsequently used the matrix-element 
generator MadGraph 5 [39] to compute all sgluon partial decay widths. We then estimate 
the total width and the different branching ratios into two gluons, into an associated pair 
of a top quark and a light quark and into two top quarks. The results are shown in Table 
2 for the two classes of considered scenarios. 

The branching of a light sgluon of few hundreds of GeV into a top-antitop pair is, as 
expected, kinematically suppressed compared to the other decay channels. This also holds 
for scenarios of class II where the sgluon field most of the time decays into a pair of gluons 
when it is light. For both classes of scenarios, the branching ratio of the tt decay increases 
with the sgluon mass (if kinematically allowed). However, the contributions of the dijet 
channel to the total width also become more important so that the branching into a top- 
antitop pair peaks for m a ~ 800 GeV and m a ~ 600 GeV for scenarios of class I and II, 
respectively, and then decreases for heavier sgluons. 

We also show in Table 2 the leading-order (LO) sgluon production cross sections as 
computed by the MadGraph 5 program for the LHC collider running at a center-of-mass 
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Scenario 


m a [GeV] 


r CT [GeV] 


BR(tt) 


BR(tj/ij) 




otot [fb] 


-^NLO 


T 
1 


zuu 


n 1 1 q 




oU/o 


^U/o 


yo < ou 


1 fi 
1.0 


I 


300 


1.05 




92.3% 


7.7% 


9902 


1.6 


T 
1 

II 


400 


z.ly 
0.286 


A A 0/ 

4.4 /o 
33.3% 


of? no/ 
ob.y /o 


O. / 70 

66.7% 


1640 


1.7 


T 
1 

II 


500 


O.bU 

0.715 


y.o /o 
47.8% 


70 ^07 

( y.o/o 


in 10/ 

lO.l/o 
52.2% 


368.5 


1.8 


T 
1 

II 


600 


A K 

4.50 

1.24 


12 /o 
48% 


(0 70 


lo/o 
52% 


97.9 


1.8 


T 
1 

II 


700 


P. OQ 

D.ZO 

1.85 


16. 1 /o 
44.7% 


7n c;oy 


lD.O/0 

55.3% 


30.2 


1.9 


T 
1 

II 


800 


7 70, 
1 . til 

2.52 


lo.O /o 
41% 


oo.y /o 


1 Q fi /, 

i y . d /o 

59% 


9.80 


2.0 


I 

II 


900 


9.43 
3.45 


13.5 % 
36.9% 


63.4% 


23.1% 
63.1% 


3.58 


2.1 


I 

II 


1000 


11.2 
4.47 


13.2 % 
33.2% 


60.2% 


26.6% 
66.8% 


1.29 


2.2 



Table 2. Dependence on the sgluon mass m a of the sgluon total width (T a ), of the different 
branching ratios (BR) to Standard Model colored particles and of the total cross section at leading 
order (otot)- The next-to-leading order if -factors i^NLO are also indicated, extracted from Ref. 
[42] . In the notations of the table, the symbol j stands for an up or charm quark while g refers to 
a gluon. 



energy of 8 TeV. The presented results correspond to the convolution of the tree-level matrix 
elements related to the Feynman diagrams of Figure 1 with the LO set of the CTEQ6 parton 
density fit [40] and renormalization and factorization scales fixed to the transverse mass of 
the produced heavy particles. As sgluon-pair production cross sections are known at the 
next-to-leading order (NLO) accuracy within the MadGolem setup [42], we also indicate 
in Table 2 the corresponding NLO if-factors. In our phenomenological analyses, signal 
event samples are then normalized according to the NLO results. 

3 Technical setup for the Monte Carlo simulations 
3.1 Objects definitions and detector simulation 

Detector simulation is performed with the DELPHES program, using the public ATLAS card 
[43]. Jets are hence reconstructed by means of an anti-Zcf algorithm with a radius parameter 
set to R = 0.4, as provided by the FAST Jet package [44]. However, jets defined with such 
a low radius, and in particular those with a small transverse momentum, may suffer from 
large bias in energy reconstruction due to the magnetic field as simulated by DELPHES, 
which spreads out the energy within the detector. Such an effect may be described by the 
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Figure 1. Tree-level Feynman diagrams associated with sgluon pair production at hadron colliders. 
These diagrams correspond to the interactions included in the Lagrangians of Eq. (2.1) and Eq. 
(2.2) and have been created by means of the program FeynArts [41]. 



variable 

^n(reco) (truth) 

"= t ^ — > 

Illrp 

where E^ ec °^ is the reconstructed jet transverse energy and E^ luth ^ is the true transverse 
jet energy denned as the jet transverse energy before detector simulation. The evolution 
of the w-variable with the (true) jet energy is presented on Figure 2 (red squares) in the 
context of dijet events issued from the decay of a sequential Z'-boson whose mass varies in 
the range [200, 1000] GeV. The energy loss reaches about 5% for jets with E^ ruth ^ = 20 GeV 
while it stabilizes at about 1% for jets with a transverse energy 

4 trUth) > 500 GeV. 

To account for this effect, an ad-hoc calibration is estimated from the above-mentioned 
dijet events. We apply on the reconstructed jet energy the correction function 

(3.2) 



^M) _ _E^ rcco ) x 



2.62. io-3 - a4Mfev m:ggl 



where is the jet transverse energy after calibration and we show the associated effects 

on Figure 2 (blue circles). This calibration procedure allows us to recover the correct jet 
energy for transverse energy as low as E^ Tuth ^ of about 40 GeV. 

In our analysis, only jets with a transverse energy Ej- > 20 GeV (after calibration) 
and a pseudorapidity \r/\ < 2.5 are retained. In addition, we estimate a 6-tagging efficiency 
of about 60%, together with a charm and light flavor mistagging rate of 10% and 1%, 
respectively. 

Charged leptons candidates 1 are requested to have a transverse momentum px larger 
than 20 GeV and a pseudorapidity \r/\ < 2.47 and \rj\ < 2.5 for electrons and muons, 



: By the generic terminology charged leptons, we only consider electrons and muons. 
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Figure 2. Evolution of the w-variable denned in Eq. (3.1) with respect to the true transverse 
energy of the reconstructed jet before (red squares) and after (blue circles) calibration. 



respectively. We also impose two isolation criteria. First, the closest jet to an electron is 
removed if the angular distance AR = \J A(p 2 + At] 2 < 0.1, where </> stands for the azimutal 
angle with respect to the bean direction. Secondly, in the case at least one jet is present 
within a cone of radius R = 0.4 centered on the lepton, the lepton is removed from the 
event . 

3.2 Monte Carlo simulations of the signal and background processes 

We perform in this work a prospective phenomenological analysis of multitop events induced 
by the production of a pair of sgluon fields. We aim to estimate the LHC sensitivity to 
the presence of sgluons in such events by means of a phenomenological analysis employing 
Monte Carlo simulations. Our analysis context consists of the LHC collider running at a 
center-of-mass energy of y/s = 8 TeV and we normalize our event samples to an integrated 
luminosity of 20 fb _1 , which corresponds to the expectations for the end of the 2012 LHC 
run [45]. 

The hard scattering processes related to both the signal and the different sources of 
background have been described with the matrix-element generator MadGraph 5 [39]. 
Using the QCD factorization theorem, the matrix elements are convolved with the leading 
order set of the CTEQ6 parton density fit [40] and the renormalization and factorization 
scales are fixed to the transverse mass of the produced heavy particles. Parton- level events 
are then integrated into a full hadronic environment by matching the hard scattering matrix 
elements to the full parton showering and hadronization infrastructure as provided by the 
PYTHIA 6 package [46]. We have then included a fast detector simulation as performed 
by the program DELPHES 1.9 [43], using, as stated above, the publicly available ATLAS 
detector card. 

The procedure above is however known not to provide an accurate description of the 
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kinematical properties of jets due to an underestimation of the hard emission effects by 
parton showering algorithms which only correctly model jets in the soft and collinear limit. 
In contrast, matrix elements properly describe additional hard and widely separated radia- 
tion, but are known to break down in the soft and collinear kinematical regions. Therefore, 
we allow for the matrix elements related to our background samples to contain zero, one, 
two, etc, additional hard jets. The resulting events are then merged following the Mangano 
(MLM) procedure [47] as implemented in the MadEvent program [48]. 

In this setup, two parton-level selection criteria are imposed. Firstly, the squared jet 
measure kj> = iaia(j^> i ,p^j)Rij related to two final state partons is asked to be larger 
than a process-dependent value (/c™ m ) 2 , where pTi and pTj are the transverse momenta 
of the two partons and Rij their angular distance in the (77, <p) plane. Secondly, the jet 
measure related to a final-state and an initial-state parton is defined as the transverse 
momentum of the final-state parton and is asked to be larger than a process-dependent 
value p™ in . The events are then passed to PYTHIA and jets are reconstructed using the 
Fast Jet program [44], making use of a fc^-jet algorithm with a cut-off scale Q m . An event 
is kept only if, for each of the reconstructed jets in the event, the jet measure between the 
jet and the parton which it is originating from is smaller than Q m . In order to maintain 
the full inclusiveness of the merged sample, extra jets with respect to the original number 
of partons are allowed in the sample with the highest multiplicity. The merging scale Q m 
is process-dependent and chosen in such a way that the differential jet rate distributions 
of the merged samples are smooth. The values chosen for the parton-level selections, the 
maximum number of included hard emissions and the merging scale are given in Table 3 for 
the various background processes. Concerning the signal, no merging with matrix elements 
of a higher multiplicity has been performed as the lowest order subprocess already contains 
many hard jets. 

Focusing first on the signal, we generate, for each scenario and each sgluon mass (see 
Table 1 and Table 2), three event samples according to the top multiplicity of the final 
state, the latter being taken equal to two, three and four top quarks, respectively. Our 
event sample normalization includes the NLO A'-factor presented in Table 2. Turning 
to the background, we are planning to require, in our analysis, the presence of at least 
one isolated charged lepton. We therefore only simulate the Standard Model background 
processes presented in Table 3, which also contains, in addition to the generation parameters 
above-mentioned, the cross sections employed for the normalization of the samples and the 
associated numbers of generated events. 

We first address the simulation of weak gauge boson production in association with 
jets. We have merged event samples containing up to four additional hard jets and the 
gauge bosons are enforced to decay leptonically. Parton-level events having also been al- 
lowed to contain tau leptons, we make use of the TAUOLA program [49] in order to handle 
their decays. We have reweighted the events according to the next-to- next-to-leading order 
(NNLO) cross sections as computed by the Fewz package [50-52] with the recent set of 
parton densities CT10 provided by the CTEQ collaboration [53]. Virtual photon contri- 
butions are included where relevant, together with a parton-level selection based on the 
dilepton invariant mass mgi > 50 GeV. 
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Table 3. We present the simulated background processes, together with the set of applied parton- 
level selection criteria (£;™ m and p™ m ), the number of allowed extra hard emissions at the matrix- 
element level (n) and the matching scale (Q m ). The numerical values employed for the cross sections 
(a) are also given, together with the number of generated events (A^). For each of the background 
process, the final state contains at least one lepton £, where I stands equivalently for electrons, 
muons, leptonic and hadronic taus. The notations incl. indicates that the associated samples are 
inclusive in the decays. We refer to the text for a more detailed description and for information on 
the adopted values for the cross sections. 

We now turn to background events related to top quark production. We generate two 
distinct ti samples, one associated to the semileptonic decay of the ti pair and one related 
to its dileptonic decay. Our merging procedure includes matrix elements containing up to 
two additional jets and we have reweighted all the events according to the production cross 
section at the NLO accuracy, including genuine NNLO contributions, as predicted by the 
HATHOR program [54, 55]. Single top event generation has been split into the generation 
of three different inclusive samples, following the usual parton-level distinction between s- 
channel diagrams where the top quark is produced in association with a b quark, t-channel 
diagrams where the top quark is produced in association with a light jet, and tW diagrams 
describing the associated production of a top quark and a VT-boson. In order to maintain 
this distinction non- ambiguous, the MLM merging procedure has not been applied. This 
avoids a possible double counting over the three channels since specific diagrams with extra 
radiation could in principle belong to several of the categories, even though the kinematical 
regimes are different. The events are then reweighted according to NLO cross sections 
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including genuine NNLO contributions [56-59]. 

Concerning diboson production, we have performed a merging of matrix elements in- 
cluding up to two additional hard jets and have normalized the cross sections to the NLO 
accuracy as provided by the MCFM package [60, 61]. In our setup, we have included virtual 
photon contributions where relevant and consequently imposed a selection on the dilepton 
invariant-mass of m« > 50 GeV. 

Rare Standard Model processes where a top-antitop pair is produced in association 
with one or two gauge bosons are also considered. We generate inclusive event samples 
for the tiW, ttZ and ttWW processes and perform MLM- merging of matrix elements 
containing up to two extra jets. We normalize the samples to the NLO cross sections as 
predicted by MCFM [62] for the ttW and ttZ processes and to the LO results as returned 
by MadGraph for the tiWW process. Finally, four-top production is also simulated but 
no merging is performed since the lowest order matrix-element already contains many hard 
jets. This event sample is normalized to the LO accuracy. 

4 Probing sgluons via multitop events at the LHC 

In this paper, we wish to estimate the LHC sensitivity to the search for multitop events 
originating from sgluon decays. According to the two classes of scenarios under considera- 
tion (see Table 1 and Table 2), sgluon pair production and decay lead to three topologies 
comprised of two top quarks and two light jets (tjtj), three top quarks and one light jet 
(tjtt) and four top quarks (tttt), the symbol t denoting equivalently a top and an antitop 
quark and the symbol j a light jet or 6-jet. The first two signatures are only considered 
assuming scenarios of class I, while four-top final states are produced in the context of both 
classes of scenarios. The corresponding final states are thus characterized by a large number 
of hard jets (between four and twelve) with an important heavy-flavor content arising from 
top decays. The large branching ratio associated with top quark hadronic decays would 
naively encourage us to search for sgluon pairs in fully hadronic final states. This analysis 
would however suffer from an overwhelming multijet background whose a correct estima- 
tion requires data-driven methods. To optimize the sensitivity of our search, we therefore 
restrict ourselves to leptonic final states. We design two analyses, one dedicated to events 
containing exactly one single lepton and another one to events with at least two leptons. 

4.1 Event selection for a multilepton plus jets signature 

Events are preselected with the requirement that they contain exactly two (for the tjtj 
topology) or at least two charged leptons (for the tjtt and tttt topologies) with a transverse 
momentum > 20 GeV. The invariant mass of the pair comprised of the two leading 
leptons m£t is also requested to be larger than 50 GeV in order to be compatible with 
the parton-level selection criterion of Section 3.2. At this stage, the total Standard Model 
background contains a large part of Drell-Yan lepton pair events (98.7% and 98.2% for the 
tjtj and tjtt j tttt topologies, respectively). To reduce this background, a selection on the 
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Number of jets Number of jets 



Figure 3. Jet multiplicity distribution after selecting events with exactly (left panel) or at least 
(right panel) two leptons, missing transverse energy f T > 40 GeV and a dilepton invariant-mass 
mu > 50 GeV. Contributions arising from top-antitop (red) and Drell-Yan production (orange) are 
factorized from the rest of the background (blue) and signal distributions for the tjtj signature 
(left panel) and for both the tjtt and tttt channels (right panel) are indicated by plain and dashed 
curves. 



missing transverse energy denned as 



visible particles 



(4.1) 



is applied. By requiring > 40 GeV, we take advantage of the fact that Drell-Yan events 
are characterized by a lack of missing energy whereas the neutrinos arising from leptonic 
top decays ensure signal events to contain a sensible quantity of missing energy. 

Jets present in sgluon-pair events originate mainly from the hadronization of 6-quarks 
and light quarks issued from top decays. In contrast, the hadronic activity of background 
events consists dominantly of initial-state radiation that leads to a lower jet multiplicity 
(see Figure 3). To maintain simultaneously a good sensitivity to the sgluon signal while 
discarding a substantial part of the background, the presence of at least three, four and 
five jets with pj, > 25 GeV is demanded for the tjtj, tjtt and tttt final states, respectively. 
In addition, we benefit from the presence of jets arising from the fragmentation of long- 
lived 6-quarks in sgluon-induced events, requiring a minimum number of one, two and three 
6-tagged jets for the tjtj, tjtt and tttt topologies, respectively. 

After applying the above-mentioned selections to the preselected dileptonic events, the 
selection efficiencies for the signal are found to range from 15% to 50% for a sgluon mass 
of m a = 400 GeV and from 25% to 60% for m a = 800 GeV for the three search channels 
(derived from Table 4). The expected Standard Model background, which has been divided 
by a factor of about 400, 3000 and 100000 for the tjtj, tjtt and tttt search strategies, 
respectively, is now largely comprised of Drell-Yan and top-antitop (plus possibly one or 
two additional gauge bosons) events. To further reduce it, we apply a specific selection 
on the fraction of dileptonic events and only retain those where the leptons have the same 
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Selections 


ml = 400 GeV 


tjtj channel 
m a = 800 GeV 


Backgrounds 




= 2 with p e T > 20 GeV 


(1.27 ±0.02) 


10 3 


r 1 r I n q-i 

o.lo ± 0.31 


(1.721 ±0.002) 


10 7 




mu > 50 GeV 


(1.16 ±0.02) 


10 3 


4.75 ± 0.30 


(1.716 ±0.002) 


10 7 




$ T > 40 GeV 


(9.46 ± 0.20) 


10 2 


4.27 ±0.28 


(1.549 ±0.004) 


10 5 


N j 


> 3 with p ] T > 25 GeV 


(9.26 ± 0.20) 


10 2 


4.27 ±0.28 


(5.693 ± 0.020) 


10 4 




N b > 1 


(6.10 ±0.16) 


10 2 


2.96 ±0.24 


(4.089 ±0.011) 


10 4 




Same sign dilepton 


(2.83 ±0.11) 


10 2 


1.12 ±0.15 


(4.191 ±0.035) 


10 3 



Selections 


mi = 400 GeV 


tjtt channel 
m T a = 800 GeV 


Backgrounds 


N e > 2 with p e T > 20 GeV 


(2.92 ±0.10)-10 2 


4.98 ±0.18 


(1.722 ±0.002) 


10 7 


m u > 50 GeV 


(2.66±0.10)-10 2 


4.70 ±0.18 


(1.717 ±0.002) 


10 7 


$ T > 40 GeV 


(2.20±0.10)-10 2 


4.36 ±0.18 


(1.598 ±0.004) 


10 5 


Nj > 4 with p{, > 25 GeV 


(1.98 ±0.08) TO 2 


4.26 ±0.16 


(2.375 ±0.012) 


10 4 


N b > 2 


84.2 ±5.8 


2.00 ±0.12 


(5.950 ± 0.040) 


10 3 


Same sign dilepton 


36.8 ±3.8 


0.82 ± 0.08 


(2.860 ± 0.080) 


10 2 



Selections 


mi = 400 GeV 


tttt channel 
m T J = 800 GeV 


Backgrounds 


N t > 2 with p l T > 20 GeV 


11.44 ± 0.34 


8.36 ± 0.24 


(1.722 ±0.002) -10 7 


m u > 50 GeV 


10.51 ± 0.32 


8.00 ± 0.23 


(1.717 ±0.002)-10 7 


$ T > 40 GeV 


8.86 ± 0.30 


7.44 ± 0.22 


(1.598 ±0.004) -10 5 


Nj > 5 with p{, > 25 GeV 


7.57 ± 0.27 


6.98 ± 0.22 


(8.11 ±0.06)-10 3 


N b > 3 


1.62 ± 0.13 


2.04 ± 0.12 


(1.88 ±0.06)-10 2 


Same sign dilepton 


0.70 ± 0.08 


0.87 ± 0.08 


10.3 ± 1.5 



Table 4. Flow-charts related to the selection strategy for the tjtj (upper panel), tjtt (middle panel) 
and tttt (lower panel) topologies in the case of a multilepton analysis. We present the remaining 
number of events, together with their associated statistical uncertainties, after each of the selection 
criteria in the context of the LHC collider running at a center-of-mass energy of t/s = 8 TeV and 
for an integrated luminosity of 20 fb _1 . For the tjtj and tjtt channels, signal events are given for 
scenarios of class I and for a sgluon mass of mi = 400 (800) GeV in the second (third) column 
of the tables. Concerning the tttt channel, we present in the second (third) column of the table 
the evolution of the number of signal events for a scenario of class I (II) with a sgluon mass of 
mi = 400 GeV (mi 1 = 800 GeV) after applying each of the selections. The associated numbers of 
background events are shown in the fourth column of the tables. For a detailed description of each 
of the selections, we refer to the text. 

electric charge. The signal efficiency of such a criterion is of about 50% while only 10% and 
20% of the background events survive in the context of the tjtj and tjtt/tttt topologies. 

The numbers of events selected at each step of the analysis are indicated in Table 4 for 
two representative benchmark scenarios with sgluon masses of 400 GeV and 800 GeV, as 
well as for the background. It turns out that after all selection criteria, the background is 
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Figure 4. Distribution of the Ht variable defined in Eq. (4.2) after applying the selection strategy 
associated with the multilepton analysis presented in the text for the tjtj (left) and tjtt (right) 
topologies. We distinguish the dominant source of background associated with the production of ti 
pairs together with jets (red) from the other contributions (blue). We superimpose the correspond- 
ing curves for two signal scenarios of class I with respective sgluon masses of 400 GeV (plain black 
curve) and 800 GeV (dashed black curve). 



dominated by top-antitop events for the tjtj and tjtt search channels and by events related 
to the associated production of a top-antitop pair with one or two additional gauge bosons 
for the tttt signature. 

In our simulation setup, the multijet background, jets faking leptons and charge misiden- 
tification have not been accounted for. However, on the basis of the analysis of Ref. [63] 
where same sign dilepton events are investigated after selection criteria similar to those 
applied in this paper, these sources of background have been found to contribute in a sig- 
nificant way as the associated number of surviving events after applying all the selections 
is ten times larger. We therefore adopt a conservative approach and derive, in Section 4.3, 
two limits on sgluon-induced new physics. The first one is extracted after omitting the 
non-simulated backgrounds and the second one is obtained after multiplying the number of 
background events by a factor of ten. 

We then consider the Ht variable defined by 

Ht= Yl \\pt\\, (4.2) 

jets, leptons, missing energy 

since signal events are expected to contain more jets and leptons than background events. 
Omitting the tttt channel as its statistical significance is very poor (see Table 4 for illustra- 
tive benchmark scenarios), we present Ht distributions for the tjtj and tjtt topologies on 
the left and right panels of Figure 4, respectively. On both panels, we show curves associ- 
ated with signal scenarios of class I where the sgluon mass is set to 400 GeV and 800 GeV. 
The distributions present a steep rise once the production threshold is reached, followed by 
a large peak centered around twice the sgluon mass. We then compare these shapes to the 
corresponding background distributions, after factorizing out the dominant ti contribution 
(in red) from the rest of the background events (in blue). This suggests to probe the LHC 
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Selections 


tjtj channel 

771,,. — 4UU <^rev tti^ — oUU Kjev rjacKgiounos 


N e = 1 with p e T > 25 GeV 
$ T > 40 GeV 
M]f > 25 GeV 
Nj > 6 with p> T > 25 GeV 
N b > 1 


(1.07 ± 0.01) TO 4 48.4 ± 1.0 (2.376 ±0.003) -10 s 
(7.72 ±0.06)- 10 3 40.1 ±0.9 (6.836 ± 0.002) TO 7 
(6.49 ±0.05) TO 3 32.5 ±0.8 (6.722 ± 0.002) TO 7 
(3.91 ±0.04)- 10 3 26.4 ±0.7 (8.634 ± 0.024) TO 4 
(2.94 ±0.04) TO 3 20.4 ±0.6 (6.407 ± 0.014) -10 4 



Selections 


tjtt channel 

mi = 400 GeV mi = 800 GeV Backgrounds 


N e = 1 with p e T > 25 GeV 
$ T > 40 GeV 
Mjf > 25 GeV 
Nj > 7 with V j T > 25 GeV 
N b >2 


(6.11 ±0.11)-10 2 11.25 ±0.20 (2.376 ±0.001) -10 s 
(4.45 ±0.09)- 10 2 9.59 ±0.18 (6.836 ± 0.002) -10 7 
(3.87 ± 0.09) -10 2 8.01 ±0.16 (6.722 ± 0.002) -10 7 
(2.04 ± 0.06) -10 2 5.86 ±0.14 (2.613 ± 0.012)-10 4 
(1.00 ± 0.05) -10 2 3.17 ±0.10 (9.330 ± 0.050) -10 3 



Selections 


tttt channel 

mi = 400 GeV m T J = 800 GeV Backgrounds 


N e = 1 with p e T > 25 GeV 
$ T > 40 GeV 
M]f > 25 GeV 

Nj > 8 with p? T > 25 GeV 
iV b >2 


34.9 ± 0.6 24.6 ± 0.4 (2.376 ± 0.001) • 10 8 
27.5 ± 0.5 21.3 ± 0.4 (6.836 ± 0.002) • 10 7 

23.8 ±0.5 18.1 ±0.4 (6.722 ± 0.002) -10 7 

10.9 ± 0.3 13.0 ± 0.3 (7.020 ± 0.060) • 10 3 
7.28 ± 0.27 8.96 ± 0.24 (2.658 ± 0.026) • 10 3 



Table 5. Same as Table 4 but in the context of a single lepton final state. For a detailed description 
of each of the selection criteria, we refer to the text. 



sensitivity to the presence of sgluon fields coupling dominantly to top quarks by means of a 
careful investigation of the shape of the entire Ht distributions, keeping all events, rather 
than requiring (unefficient) selection on this variable. Limits at the 95% confidence level 
are extracted in this way in Section 4.3. 

4.2 Event selection for a single lepton plus jets signature 

Final states containing exactly one lepton are expected to be copiously produced from 
sgluon pair production and decay at the LHC as they correspond to a fraction of events 
equal to 36% and ~ 41% in the context of the tjtj and tjtt /tttt topologies, respectively. We 
preselect events by requiring exactly one single lepton with a transverse momentum p l T > 25 
GeV. Consequently, the Standard Model background is expected to be dominated at 92% by 
events associated with the production of a VF-boson in association with jets. The expected 
number of sgluon events ranges from 34.9 (24.6) to 10700 (48.4) for a representative scenario 
with a sgluon mass of 400 GeV (800 GeV), as shown in Table 5. 

The multijet background has not been taken into account up to now. However, data- 
driven methods used to estimate its shape and normalization tend to show that it may 
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Number of jets 



Figure 5. Jet multiplicity distribution after selecting events with exactly one lepton, missing 
transverse energy $ T > 40 GeV and a W-boson transverse mass Mjf > 25 GeV. We distinguish 
the tt (red) and W-boson plus jets (orange) contributions from the rest of the background (blue) 
and present signal distributions for a sgluon scenarios of class I with m a = 400 GeV in the tjtj 
(plain), tjtt (strong dashed) and tttt (light dashed) channels. 



have to be considered in our analysis [64]. In order to realistically and reliably reject this 
background, we impose, inspired by the experimental analysis of Ref. [64], the missing 
transverse energy of the events to be larger than $t > 40 GeV and the reconstructed 
VF-boson transverse mass, defined as 



= ^2jfi T $ T [l - cos A^ T j , (4.3) 

to be larger than 25 GeV. In this equation, the quantity Acj)^ g stands for the angular 
distance, in the azimutal direction with respect to the beam, between the lepton and the 
missing energy. 

As for the multilepton analysis, a higher multiplicity of hard jets is expected in signal 
events. They indeed arise from the hadronization of final state quarks, in contrast to the W- 
boson plus jets background where jets originate predominantly from initial-state radiation 
(see Figure 5). Events are therefore selected with the requirement that they contain at 
least six, seven and eight jets with a transverse-momentum p J T > 25 GeV for the tjtj, tjtt 
and tttt topologies, respectively. Moreover, sgluon events are expected to include a higher 
number of heavy-flavored jets. The minimal number of 6-tagged jets is therefore required to 
be larger that one and two for the tjtj and tjtt /tttt search channels, respectively. At this 
stage, the expected Standard Model background is composed mainly of ti events, where the 
top-antitop pair is possibly produced in association with one or several gauge bosons. 

The number of events surviving to each of the selection criteria is presented in Table 
5 for two representative signal scenarios with m a = 400 GeV and 800 GeV as well as for 
the background. After all selections, we expect a number of signal events ranging from 7.28 
(8.96) to 2940 (20.4) for a sgluon mass of 400 (800) GeV according to the search channel 
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Figure 6. Distribution of the Ht variable defined in Eq. (4.2) after the selection strategy associated 
to a single lepton analysis as presented in the text. We distinguish the background events associated 
to the production of tt pairs together with jets (red) from the other contributions (blue) . For the 
tjtt channel (left panel), we superimpose the corresponding curves for two signal scenarios of class 
I with respective sgluon masses of 400 GeV (plain black curve) and 800 GeV (dashed black curve). 
In the case of the tttt search-channel (right panel), we rather consider a scenario of class I with a 
sgluon mass of 400 GeV (plain black curve) and a scenario of class II with a sgluon mass of 800 
GeV (dashed black curve). 



under consideration. In contrast, the Standard Model predicts a background of 64070, 9330 
and 2658 events for the tjtj, tjtt and tttt topologies. 

Due to the larger activity expected in signal events compared to background events, we 
consider, as in Section 4.1, the Ht variable as a discriminant between signal and background. 
Following the definition of Eq. (4.2), we present, in Figure 6, the Ht distributions after 
following the tjtt (left panel) and tttt (right panel) single lepton selection strategy intro- 
duced above. We show background distributions after distinguishing events associated with 
top-antitop production in association with jets (red) from the rest of the background (blue). 
The resulting behavior underlines a steep rise once the top-antitop production threshold 
is reached, both in the tjtt and tttt channels, followed by a peak around Ht ~ 500 GeV 
and a smooth fall with increasing energy. In comparison, signal distributions are associated 
with a clear peaky behavior centered around an Ht value depending on the sgluon mass 
(~ l.Smo-) and the tail of the distributions does not extend to a very large hadronic energy, 
in contrast to the Standard Model background expectation. Instead of requiring a selection 
criterion based on the Ht variable, we then employ the shape of the associated distributions 
as inputs to extract the LHC sgluon mass reach in Section 4. 

In the case of a tjtj signature with exactly one final state lepton, the reconstruction 
of the sgluon mass is possible if the missing energy is assumed to originate only from the 
neutrino issued from the VF-boson decay. Assigning the labeling of the six jets according 
to the process 

pp -> aa -> (tj 5 )(tj 6 ) ->■ (jii2j3j5)(i4^^j6) , (4.4) 
all possible permutations of the six jets are performed and the one minimizing a x 2 -variable 
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is retained as the true configuration of a given event. The analytical expression of this 
% 2 -variable exactly mimics Eq. (4.4). The three jets j\, j'2 and j'3 are the ones issued from 
the hadronically decaying top quark, j'3 being a 6-jet 2 , and this information is encompassed 
into the first two terms of the x 2 variable. The invariant mass of the two light jets ji and 
j2, denoted by rrij 1 j 2 , is asked to be compatible with the VF-boson mass and the three-jet 
invariant mass rrij 1 j 2 j 3 is required to be compatible with the top mass. However, in order not 
to introduce a correlation between the first two terms of the x 2 variable, we subtract from 
the reconstructed top mass TUj 1 j 2 j 3 the reconstructed dijet invariant-mass rrij 1 j 2 . The values 



of the fit parameters are taken as = 80.7 GeV, OyJ = 8.9 GeV, = 90.8 GeV and 

(r) 

o\yj = 10.5 GeV. These numerical values have been extracted from a fit based on the Monte 
Carlo truth, which ensures that each reconstructed object is correctly assigned according 
to the configuration of Eq. (4.4). In our simulation setup, it must be noted that the large 
values of 0(10%) for the widths are mainly related to detector resolution. 

In the third term of Eq. (4.5), we focus on the leptonically decaying top quark and 
ask the invariant mass rri£ u j 4 to be compatible with the top mass. Since the neutrino four- 
momentum is entirely reconstructed from the assumption that both the missing energy and 
the charged lepton are issued from the decay of a VF-boson, there is no need to add an extra 
term in the definition of the x 2 , the information being already implicitly included in the 
mevjt term. From the Monte Carlo truth, we have calculated a reconstructed mass equal 



(r) 



,(r) 



to m 



(r) 

tl 



167.8 GeV, together with a width of = 19.1 GeV. 



Finally, the definition of the x 2 variable of Eq. (4.5) exploits the fact that both the tj's 
and tj% pairs are issued from the decay of a sgluon field. The corresponding invariant masses 
5 and m£ u j i j 6 must hence be self-compatible, up to the detector resolution. This is 



313233,35 



translated in the last term of Eq. (4.5), after subtracting the respective reconstructed top 
masses to avoid possible correlations among the different terms of the x 2 . From the Monte 



Carlo truth, we fix the related width to a 



(r) 
at 



0.098 GeV. 



Performing an event selection as described in the upper panel of Table 5 and minimizing 
the x 2 variable of Eq. (4.5), each event is reconstructed according to the pattern given in 
Eq. (4.4). This allows us to extract the sgluon mass ■ The resulting distributions are 
presented on Figure 7 for both the background and four signal scenarios of class I with 
a sgluon mass parameter m a taken equal to 200 GeV, 400 GeV, 600 GeV and 800 GeV, 
respectively. As expected, the background distribution, where we again distinguish the top- 
antitop pair contributions (in red) from the other sources of background (in blue), presents 
a rising behavior once the production threshold of a top-antitop pair is reached, quickly 



We are however not making use of jet flavor information in our kinematical fit. 
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Figure 7. Reconstructed sgluon mass rtv a for both the background and several signal scenarios in 
the context of the tjtj topology. Event selection is performed as in the upper panel of Table 5 and we 
reconstruct each event according to the pattern given in Eq. (4.4) by means of the minimization of 
the x 2 variable of Eq. (4.5). Concerning the background distribution, we distinguish the dominant 
top-antitop contribution (red) from the rest of the background (blue) and superimpose signal curves 
for four scenarios of class I, with respective sgluon mass of 200 GeV (plain) , 400 GeV (strong dotted) , 
600 GeV (light dotted) and 800 GeV (dash-dotted). 



followed by a slow fall which extends to rather large values of the reconstructed sgluon mass 
(r) 

m& . In contrast, the signal distributions all show a peak. For light sgluons (rn a = 200 
GeV or 400 GeV), this peak is clearly centered around the true sgluon mass. In the case 
of heavier sgluons (rn a = 600 GeV and 800 GeV), the detector resolution renders the peak 
very wide and centered around a value equal to about 70% — 80% of the real sgluon mass. 

(r) 

We then use, in the next section, the reconstructed mass m a to extract limits on the sgluon 
mass reachable at the LHC. 

4.3 LHC sensitivity to a sgluon field dominantly coupling to top quarks 

For each of the final states considered in this paper, we combine the number of expected 
signal and background events, including their corresponding statistical uncertainties, to 
calculate upper limits on the signal cross section at the 95% confidence level. To this aim, 
we use the CLs technique [65] as implemented in the MCLlMIT software [66]. We employ the 
Ht variable defined in Eq. (4.2) to discriminate signal from background in the multilepton 
(for the tjtj, tjtt and tttt topologies) and single lepton (for the tjtt and tttt topologies) 
analyses while the reconstructed sgluon mass is chosen in the case of a single lepton analysis 
applied to a tjtj final state 3 . 

3 We have checked that considering the reconstructed mass instead of the Ht variable allows us to 
improve the LHC sensitivity by about 15%-20% in the low mass region. This choice is however irrelevant 
for the higher sgluon mass region. 
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Single lepton analysis 


Multilepton analysis 


Multilepton analysis 
(background xlO) 


tjtj 
tjtt 

tttt (Scenario I) 
tttt (Scenario II) 


590 tf GeV 
480 tlo GeV 

640 tto GeV 


570 tf GeV 
520 jjjj GeV 

650 tlo GeV 


440 ±f° GeV 
520 tno GeV 



Table 6. Expected sensitivity of the LHC collider, running at a center-of-mass of 8 TeV and for 
an integrated luminosity of 20 fb _1 , to sgluons. The results are given, together with the associated 
lcr statistical uncertainties, in terms of upper bounds on the sgluon mass to be reached, at the 
95% confidence level, for different types of analyses and scenarios. For scenarios of class I, the 
expectations for the tjtj channel are given on the first line of the table, those related to the tjtt 
channel on its second line and those associated with the tttt channel on its third line. Scenarios of 
class II are only investigated in the context of the tttt search-channel and results are shown on the 
fourth line of the table. Finally, we distinguish results obtained by employing a single lepton analysis 
(second column), a multilepton analysis after neglecting the QCD background (third column) and a 
multilepton analysis after multiplying the simulated background by a factor of ten (fourth column) . 

The results are presented in Figure 8 for the tjtj channel (upper panel), tjtt channel 
(middle panel) and tttt channel (lower panel) as dashed and dot-dashed curves in the context 
of the multilepton and single lepton analysis, respectively. In addition, the effects of the non- 
simulated sources of background mentioned in Section 4.1 are presented as dotted curves. 
On these figures, we also show the theoretical cross sections related to sgluon-induced 
production of multitop final states as a function of the sgluon mass. In addition to the 
central NLO curves derived from Table 2, we include a 30% uncertainty band corresponding 
to the typical order of magnitude of the variations of the results with respect to different 
choices for the factorization and renormalization scales [42]. The curves associated with 
scenarios of class I are shown in dark gray, while those related to scenarios of class II are 
given in light gray. 

Sgluon masses excluded at the 95% confidence level are also indicated in Table 6, for 
each final state and for each scenario considered in this paper. The quoted uncertainties 
correspond to an estimation of the cross section limits at the la level, where a stands for 
the statistical uncertainty. 

In the multilepton analysis, sgluon masses lower than 570 and 520 GeV can be excluded 
in the tjtj and tjtt topologies, respectively, for scenarios of class I. Equivalently, the ATLAS 
experiment is sensitive to sgluon-induced multitop production cross sections of 0(100) fb 
for both the tjtj and tjtt signatures at \fs = 8 TeV. Inspecting the second column of 
Table 6, the limits obtained for the tjtj search-channel barely vary (by less than 10%) 
when accounting for la statistical uncertainties. In contrast, the bounds extracted from 
the analysis of the tjtt topology are found to be more sensitive to statistics as they can 
vary by about 17% with respect to (un)lucky fluctuations. This feature can be understood 
from the behavior of the theoretical cross sections in the 500-600 GeV sgluon mass range. 
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Figure 8. The 95% confidence level expected signal cross sections as a function of the sgluon mass, 
for an integrated luminosity of 20 fb -1 and s/s = 8 TeV. The bands around the theoretical curves 
correspond to an uncertainty of 30%, and predictions for scenarios of class I (II) are presented as 
light (dark) gray bands. The expected limits in the tjtj (upper), tjtt (middle) and tttt channels 
(lower) are given, for the single lepton (multilepton) analysis, as dot-dashed (dashed) curves. In 
the multilepton case, we also show the results obtained after enhancing the background by a factor 
of ten (dotted curves) . 
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While the tjtj cross section decreases with the sgluon mass, the tjtt cross section is rather 
flat. The multilepton analysis is however not sensitive to sgluon-induced production of four 
top quarks for class I scenarios, at sfs = 8 TeV and for an integrated luminosity of 20 fb" 1 . 
This results from the too low branching ratio of the sgluon decay into a top-antitop pair 
(see Table 2). In class II scenarios, this branching ratio is 2.5 — 7.6 times more important, 
so that sgluon masses lower than 650 GeV, or cross section of O(10) fb, can be excluded at 
the 95% confidence level. 

When taking into account non-simulated backgrounds, the total number of expected 
background events is conservatively multiplied by ten. Consequently, the tjtt analysis is 
not sensitive to sgluon-pair production anymore, while the masses excluded when analyzing 
tjtj and tttt (for scenarios of class II) final states decrease from 570 to 440 GeV and 650 to 
520 GeV, respectively. 

We now turn to single lepton analyses and show that the reconstructed sgluon mass 
from a tjtj topology can be used to exclude, at the 95% confidence level, sgluon as heavy 
as 590 GeV (dot-dashed curve). Moreover, these bounds are found not to drastically vary 
when including la statistical uncertainties. Concerning the tjtt and tttt topologies, the Ht 
variable is considered as a discriminant and sgluon masses up to 480 GeV and 640 GeV can 
be reached. In the first case, statistical fluctuations can lead to different expectations by 
about ±15% in the tjtt case while in the tttt case, the results are found only to slightly 
change by about 5%. 

5 Conclusion 

Many new physics theories predict the existence of a scalar field, commonly dubbed sgluon, 
lying in the adjoint representation of the QCD gauge group. To investigate the sensitivity 
of the LHC to this particle in the case it couples dominantly to the top quark, as motivated 
by hybrid N =1/N = 2 or i?-symmetric supersymmetric theories, an effective field theory has 
been built. This theory consists of a minimal extension of the Standard Model including 
a single real sgluon field, together with a set of interactions leading to its production and 
decay at the LHC. The model has been implemented in the FeynRules package so that a 
UFO model for MadGraph has been extracted. 

Final states containing two, three, or four top quarks produced in association with jets 
have been investigated. To optimize the sensitivity of the search for sgluon-pair production 
predominantly decaying to top quarks, two signatures (multilepton and single lepton) have 
been considered in three topologies (tjtj, tjtt, tttt). We have shown that sgluons of about 
500-700 GeV can be reached at the LHC collider running at a center-of-mass energy of 
8 TeV and for an integrated luminosity of 20 fb" 1 . Equivalently, this mass constraint can 
be translated into a bound on the cross section associated with the production of two, 
three, or four top quarks issued from the decays of a pair of sgluons. Contributions higher 
than 10-100 fb can be excluded for most scenarios, at the 95% confidence level. While 
a search strategy based on a single lepton selection ensures, with appropriate selection 
criteria, the multijet background to be under good control, non-simulated contributions 
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to the background such as jets faking leptons and charge misidentification can reduce the 
sensitivity of the multilepton signature by several hundreds of GeV. 

Our work therefore motivates a future extension of the two performed analyses in the 
context of a full detector simulation of the LHC experiments. 
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